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ABSTRACT 


In this study Pickover biomorphs are analysed as being dependent on the chosen complex number system 
in which iterations of analytic functions are performed. Moran’s spatial autocorrelation function and two 
forms of entropy, the Shannon entropy and the sample entropy, are chosen in order to find correlations 
and measure complexity in Pickover biomorphs. These turn out to be strongly correlated and low-entropy 
objects with a fractal dimension between 1.4 and 2. It is shown that there is a strong maximum in 
Keywords: correlation and a strong minimum in entropy for the case of Galilean complex numbers corresponding to 

Pickover biomorphs the square of the generalised imaginary unit being equal to zero. 
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1. Introduction 

Mandelbrot sets, Julia sets [1-4] and Pickover biomorphs 
[5-9] belong to a group of mathematical and graphical fractal 
objects which are of particular fascination to both scientists and 
artists. Remarkably some Pickover biomorphs morphologically re¬ 
semble unicellular protozoa with their various organelles or even 
aquatic crustaceans. Together with better known related objects 
such as Mandelbrot sets, Pickover biomorphs are an interesting ex¬ 
ample of very complex structures emerging from the application of 
simple rules which determine the application of a transformation 
of the complex plane by iterations of analytic functions. Due to the 
simplicity of such rules, Pickover biomorphs have attracted some 
mild attention in the field of theoretical biology [10,11]. These ap¬ 
plications in the simulation of evolution set out the most impor¬ 
tant context of the present study. In one sense, Pickover biomorphs 
may also be said to be a fascinating by-product of complex dynam¬ 
ics in mathematics [12]. 

It is not, however, easy to analyse Pickover biomorphs mathe¬ 
matically partly because (to the authors’ best knowledge) there is 
not even a heuristic mathematical definition. Pickover biomorphs 
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are in fact defined by a family of algorithms which create them 
or even directly as merely graphical objects which resemble cer¬ 
tain invertebrates. Thus, this study solely relates to images ob¬ 
tained from more or less standard algorithms. In order to charac¬ 
terise biomorphs quantitatively, two kinds of measures have been 
used. Firstly, as a measure of spatial correlations, the Moran au¬ 
tocorrelation function has been chosen. Secondly, as a measure of 
complexity, two versions of entropy, the Shannon entropy and the 
sample entropy have been applied. The second type of entropy is 
considered to be an approximation of the Kolmogorov-Sinai en¬ 
tropy. It turns out that even though the autocorrelation function 
and the entropy measure different characteristics of the biomorphs 
they still seem to be complementary to each other. 

There are two obvious considerations at this point, namely, 

(a) whether it is of any use to extend the notion of biomorphs 
to non-standard systems of complex numbers, and (b) what ex¬ 
actly the purpose of application of the machinery of statistics and 
complexity-related measures to them is. For (a), we believe that 
from the point of view of investigations in evolution theory exem¬ 
plified by Mojica et al. [10] it is useful to have at one’s disposal 
structures with sharply defined qualitative changes of properties 
which can be obtained just by changing one very well defined pa¬ 
rameter, in our case here the square of the imaginary unit. For 

(b) , when simulating the evolution of simple forms with the help 
of genetic algorithms as in [10] it is important to know whether 
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the “evolution” leads to greater or smaller complexity, stronger or 
weaker morphological correlations, and even lower or higher frac¬ 
tal dimensions. 

2. Non-standard systems of complex numbers 

It is well known that one of the fundamentally important pro¬ 
cedures in mathematics is the pairing of two real numbers to form 
a complex number. Multiplication in the set of complex numbers 
can be defined using Hamilton’s construction or alternatively, any 
complex number z can be written as z = x + iy with i 2 = -1. One 
can then formally operate with complex numbers as with real ones 
under the assumption of distributiveness of multiplication with re¬ 
spect to addition and substituting in for i 2 . Let us note that if i is 
understood as just a symbol (rather than a number) with a defined 
square, the construction of complex numbers can easily be ex¬ 
tended to cases where i 2 is different to -1. In fact there exist three 
canonical complex number systems corresponding to i 2 = -1, 0 
and +1. Following Hucks [13], these will be referred to here as 
standard (Euclidean), Galilean, and hyperbolic (complex) numbers 
respectively. In the case of Galilean complex numbers [14] the law 
of multiplication of two numbers is: 

(a + bi)(c-\-di) = ac + (ad + bc)i, (1) 

and in the case of hyperbolic complex numbers the law takes the 
form: 

(a + bi) (c + di) = (ac + bd) + (ad + bc)i . (2) 

An important point here is that the concept of analytic (or holo- 
morphic) function as a function which depends only on z = a + bi 
but not on z* = a - bi can be naturally extended to the case of 
Galilean and hyperbolic numbers. The Cauchy-Riemann equations 
as necessary conditions of analyticity can also be easily formulated. 

It can be noted that hyperbolic complex numbers seem to 
be particularly interesting from the point of view of theoretical 
physics. Indeed, a reasonably far-reaching theory of analytic func¬ 
tions can be built based on their set. The corresponding Cauchy- 
Riemann equations are no more than the free Maxwell equa¬ 
tions in (1+1) dimensions. Taking into account that every ana¬ 
lytic function provides solutions to the Cauchy-Riemann equations, 
solutions to (1+1) -dimensional Maxwell equations can be ob¬ 
tained, which might satisfy interesting boundary conditions [13]. 
This could be of interest, for example, from the point of view of 
cavity quantum electrodynamics with moving boundaries where 
exact oscillating solutions to the Maxwell equations are important 
(see, e.g. [15-17]). 

3. Results 

In the first part of this section, several pictures of biomorphs 
obtained from the iterations of almost the simplest possible an¬ 
alytic functions, namely quadratic and cubic ones, are presented. 
The second part is devoted to their analysis using well-defined sta¬ 
tistical and complexity measures. In the latter part, the study has 
been deliberately restricted to such simple functions in order to 
avoid dealing with too much scattered data with the loss of neces¬ 
sary focus. 

3.1. Standard Pickover biomorphs from iteration of a cubic function 

Let us consider the following iteration: 

z <— z n + C, (3) 

where <— is the assignment operator (as in [18]) and we then in¬ 
troduce the following criterion to stop the iteration procedure: 

Re(z) 2 > r 0 or Im(z ) 2 > r 0 . (4) 


It is precisely because of the seemingly innocent “or” statement in 
the above criterion that the images produced by the appropriate 
colouring of the structure obtained differ radically from otherwise 
closely related Mandelbrot or Julia sets. Examples of standard Pick- 
over biomorphs obtained from the above procedure with n = 2 and 
n = 3 are shown in Fig. 1. The “simulation cell” has a dimension 
of 1024x1024 pixels and the starting points have been located in 
the square -2.5 < Re(z ) < 2.5, -2.5 < Im(z) < 2.5. In all examples 
shown below the maximal number of iterations was equal to 100. 

The value of r 0 has been set equal to 10.0 and C has been set 
equal to 0.5 (i.e. with a zero imaginary part). 

3.2. Pickover biomorphs in the realm of Galilean numbers 

The same quadratic and cubic iteration schemes and parameters 
are applied to Galilean complex numbers and result in the image 
shown in Fig. 2: 

3.3. Pickover biomorphs and hyperbolic complex numbers 

If it is now assumed that i 2 = +1, the Pickover biomorphs take 
very different forms as illustrated in Fig. 3. 

The self-similarity though evident is clearly restricted and simi¬ 
lar structures seem to be shifted and spatially separated from each 
other. 

3.4. Pickover biomorphs for the squares of imaginary units dependent 
on points in the complex plane 

The study has also experimented with a type of complex num¬ 
bers with the following multiplication law: 

(a + bi)(c + di) = (ac + bdi 2 (a, b, c, d)) + (ad + bc)i, (5) 

with the obvious requirement of symmetry i 2 (a, b, c, d) = 
i 2 (c,d,a,b ). The following images have been obtained from the 
iterations of a cubic function and three different Gaussian forms of 
i 2 depending only on a and c. On the left-hand side of this figure 
i 2 = -1 + exp(-2(a 2 + c 2 )), at the centre i 2 = exp(-2 (a 2 + c 2 )), 
and on the right-hand side i 2 = -1 + 2exp(-2(a 2 + c 2 )) (Fig. 4). 

3.4.1. Correlations in images of Pickover biomorphs 

As a measure of correlations in the images of Pickover 
biomorphs Moran’s autocorrelation function / [19 has been cho¬ 
sen and is defined as: 

. _ N EiJ w ij ( x i - x ) (Xj - x) 
w E+i-x) 2 ’ 

where N is the number of spatial units indexed by i and j, x is the 
mean of all {xj, Wy is a matrix of weights such that w ti = 0, and 
W = J2i j w ij- ^ is well known that the value of I may strongly de¬ 
pend on the choice of matrix of weights. An attempt has not been 
made to obtain any comprehensive picture of such dependence in 
our case, but instead a simple binary matrix has been applied. 
Firstly, Wy was chosen to be 1 for the four nearest neighbours 
and 0 otherwise. Secondly, the computations were then repeated 
for Wfj = 1 for the six nearest neighbours. These choices of neigh¬ 
bours correspond to the von Neumann and Moore neighbourhoods 
in cellular automata theory [20]. The biomorph for analysis was 
obtained using a grey-scale picture of the size of 1024 x 1024 pix¬ 
els. 

It appears that the correlations in biomorphs do not depend 
much on the system of weights in the Moran autocorrelation func¬ 
tion. Furthermore, the qualitative outcome of the dependence on i 2 
is the same for quadratic and cubic functions. We have observed in 
addition that Moran’s function grows when the size of the pictures 
is increased from 128 x 128 to 1024 x 1024 (Figs. 5 and 6). 
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Fig. 1 . Standard Pickover biomorphs obtained from the iteration of quadratic (left) and cubic (right) functions. 



Fig. 2. Pickover biomorphs for Galilean complex numbers obtained from the iteration of quadratic (left) and cubic (right) functions. 



Fig. 3. Pickover biomorphs for hyperbolic complex numbers. 




















A. Jakubska-Busse et al./ Chaos, Solitons and Fractals 113 (2018) 46-52 


49 



Fig. 4. Pickover biomorphs obtained from the iteration of a cubic function for i 2 depending on points in the complex plane. 




Fig. 5. Moran’s I global autocorrelation function as dependent on the square of the imaginary unit for a quadratic function. Left: Matrix of weights with 4 nearest neighbours: 
Right: Matrix of weights with 6 nearest neighbours. 




Square of imaginary unit Square of imaginary unit 


Fig. 6. Moran’s I global autocorrelation function as dependent on the square of the imaginary unit for a cubic function. Left: Matrix of weights with 4 nearest neighbours: 
Right: Matrix of weights with 6 nearest neighbours. 


3.5. Entropy as a measure of complexity in images of biomorphs 

To compute the global sample entropy [21] as well as the Shan¬ 
non entropy of a biomorph, the matrix of grey-scale representation 
of a given image was used and then flattened to obtain a one¬ 
dimensional representation of the biomorph. This has been done 
for 50 values of the square of the imaginary unit in order to obtain 
a fairly comprehensive understanding both of the generic value of 
the entropy and its dependence on i 2 . The results are shown in 
Figs. 7 and 8. 


In the calculations of sample entropy we have followed the 
standard prescriptions regarding the length m of the data segment 
being compared (taken equal to 2) and similarity criterion r (being 
equal to one fifth of the standard deviation). The sample entropy 
tends to be smaller for larger values of i 2 , apart from the obvious 
local minimum near zero. Thus, the “time series” which represent 
the biomorph appear to have more repetitive patterns for those 
larger values. It is interesting to note that the Shannon entropy 
does not distinguish between the signs of i 2 and remains almost 
flat apart from the region close to i 2 = 0. 
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square of imaginary unit 

Fig. 7. Sample entropy as a function of the square of the imaginary unit. Left: for the quadratic function; Right: for the cubic function. 



3.6. Hurst exponent 

The Hurst exponent is one of the measures of the long-term 
memory of a time series [22,23]. It can be used as an indicator 
of irregularity, especially of time-series. Processes with higher (i.e. 
closer to 1) Hurst exponent appear to be more regular and indicate 
the presence of long-time correlations, [24]. Extension of the con¬ 
cept to any sequence of values (not necessarily labelled with time) 
is straightforward. 

In our case, values of the Hurst exponent using the rescaled- 
range R/S algorithms were also obtained through the use once 
again of a one-dimensional matrix obtained by flattening the orig¬ 
inal grey-scale representation of the biomorphs. The Hurst expo¬ 
nent has been computed from the resulting “time series”. 

A visible and interesting feature of the dependence of the Hurst 
exponent on i 2 is that it is very different for two iterated func¬ 
tions, and the structure of both dependencies is rich with several 
minima and maxima. For the cubic function, for all values of i 2 
the Hurst exponent has turned out to be less than 0.5 indicating 
that the fractal dimension of the one-dimensional representation 
of the biomorph may be between 1.5 and 2 with a strong maxi¬ 


mum near i 2 = 0. But it is not true for the quadratic function, al¬ 
though the neighborhood of the point i 2 = 0 is clearly well distin¬ 
guished (Fig. 9). 

It appears that the “time series” associated with the biomorphs 
are quite close to “Brownian motion” since the Hurst exponents 
are close to 0.5. However, weak “anti-persistency” is evident in the 
case of the cubic function. We have not yet found a proper inter¬ 
pretation of the details of the behaviour of Hurst exponents since 
more cases will need to be investigated to do this. 

4. Discussion 

In this work the properties of specific, biologically-inspired im¬ 
ages obtained from iterations of complex functions have been anal¬ 
ysed. These functions have been applied to three different systems 
of complex numbers: Euclidean (i.e. standard), Galilean, and hyper¬ 
bolic complex numbers. The following conclusions can be drawn. 
Firstly, the concept of Pickover biomorphs can be extended from 
the set of complex numbers to the set of Galilean and hyperbolic 
complex numbers, and more generally to the set of numbers with 
the square of the imaginary unit being a function of points in the 















A. Jakubska-Busse et al./ Chaos, Solitons and Fractals 113 (2018) 46-52 


51 




Fig. 9. The Hurst exponent as a function of the square of the imaginary unit. A grey-scale representation of each biomorph was first obtained and then its matrix repre¬ 
sentation was flattened to obtain a one-dimensional sequence of values from which the Hurst exponent has been computed. Left: for the quadratic function; Right: for the 
cubic function. 


complex plane. Secondly, the graphical structures which emerge 
from iterations of (hyperbolic and Galilean) complex functions are 
usually simpler and less aesthetically appealing than in the case 
of standard (Euclidean) complex functions; however, self-similar 
structures are likely to appear. Thirdly, several measures of corre¬ 
lations as well as complexity such as Moran’s correlation function, 
sample entropy and the Hurst exponent can be useful in the anal¬ 
ysis of biomorphs associated with any concept of complex num¬ 
bers. Fourthly, whilst all the biomorphs are graphical objects in 
which naturally, strong spatial correlations appear, the latter are 
less pronounced in the case of standard (Euclidean) complex num¬ 
bers than in the case of hyperbolic numbers with a somewhat un¬ 
expected deep local minimum on the side of the negative values 
of the square of the imaginary unit i 2 , and a strong maximum very 
close to the zero value of i 2 . Fifthly, the dependence of sample en¬ 
tropy (obtained by flattening the matrix representation of the grey¬ 
scale image representing a biomorph) appears to be complemen¬ 
tary with respect to the Moran autocorrelation function. The sam¬ 
ple entropy is smaller for Galilean and hyperbolic numbers than 
for Euclidean ones, indicating stronger self-similarity in the former 
case. This is, however, mostly due to the higher overall symmetry 
of biomorphs obtained with non-negative i 2 . Finally, the Hurst ex¬ 
ponent has been obtained as a function of i 2 . 

5. Methods 

All programs for the purpose of this work were written in 
Python (version 2.7). The following Python modules were used: 
numpy, scipy, matplotlib, Pillow [25 , pysal [26], pyentrp [27], and 
nolds [28]. 

In explanation of the process of the generation of the figures, 
the drawing area of 1024 x 1024 pixels was chosen in order to ob¬ 
tain the drawings of the biomorphs. The parameter C which enters 
the iteration algorithm was set equal to (0.5, 0) and the drawing 
area was mapped to an area of (generalized) complex plane of the 
size ([-2.5, 2.5] x [-2.5, 2.5]). The maximal number of iterations 
has been set equal to 100; the programs left the iteration loops if 
the iteration index exceeded 100, or if either the real or the imag¬ 
inary part of z exceeded 10.0. Then the final values of the real and 
the imaginary parts of z were tested for whether any of them be¬ 
long to an interval (a, b). Depending on (a, b), a specific colour of 
the pixel corresponding to the initial value of z has been chosen as 
illustrated by the following Tables 1 and 2: 

In the above tables, “RGB” signifies the amount of the Red, 
Green, and Blue colours in the colour of a pixel. The code of 
Table 3 has also been used to produce Table 4. 


Table 1 

Colouring code for the biomorphs obtained on using standard (Euclidean) com¬ 
plex numbers (1). 


a 

b 

RGB 

0.0 

0.5 

(0, 0, 0) 

0.5 

4.0 

(0, 150, 0) 

4.0 

7 

(150, 0, 0) 

7.0 

10.0 

(0, 0, 150) 

10.0 

13.0 

(150, 150, 0) 

13.0 

17.0 

(150, 0, 150) 

17.0 

20.0 

(0, 150, 150) 

20.0 

23.0 

(150, 150, 150) 

23.0 

oo 

(255, 255, 155) 


Table 2 

Colouring code for the bio morphs obtained on using Galilean complex numbers (2). 


a 

b 

RGB 

0.0 

10-1° 

(50, 50, 50) 

10-io 

1.0 

(0, 150, 0) 

1.0 

2.5 

(150, 0, 0) 

2.5 

4.0 

(0, 0, 150) 

4.0 

5.0 

(150, 150, 0) 

5.0 

6.0 

(150, 0, 150) 

6.0 

7.0 

(0, 150, 150) 

7.0 

8.0 

(255, 0, 15) 

8.0 

oo 

(255, 255, 155) 


Table 3 

Colouring code for the biomorphs obtained 
bers (3). 

on using hyperbolic complex num- 

a 

b 

RGB 

0.0 

0.5 

(0, 0, 0) 

0.5 

4.0 

(0, 150, 0) 

4.0 

7.0 

(150, 0, 0) 

7.0 

10.0 

(0, 0, 150) 

10.0 

13.0 

(150, 150, 0) 

13.0 

17.0 

(150, 0, 150) 

17.0 

20.0 

(0, 150, 150) 

20.0 

23.0 

(150, 150, 150) 

23.0 

30.0 

(250, 250, 15) 

30.0 

oo 

(255, 255, 155) 
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